Spin-Polarized Transport through Double Quantum Dots 
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We investigate spin-polarized transport phenomena through double quantum dots coupled to fer- 
romagnetic leads in series. By means of the slave-boson mean-field approximation, we calculate the 
conductance in the Kondo regime for two different configurations of the leads: spin-polarization of 
two ferromagnetic leads is parallel or anti-parallel. It is found that transport shows some remarkable 
properties depending on the tunneling strength between two dots. These properties are explained 
in terms of the Kondo resonances in the local density of states. 

PACS numbers: 



I. INTRODUCTION 

Since the Kondo effect was observed in a quantum dot 
system[lly|, electron correlation effects in such nanoscale 
systems have attracted much attention. It was uncovered 
that the Kondo effect gives rise to remarkable proper- 
ties in electron transport, which are not expected for 
free electron systems 0, H> IM IE Q? Ill- For example, 
it makes tunneling transport possible at low tempera- 
tures even if the system is in the Coulomb blockade 
regime. Since the Kondo effect originates from many 
body effects due to internal spin degrees of freedom, 
transport properties are quite sensitive to applied mag- 
netic fields 0, Q- This fact has naturally stimulated 
intensive works on spin-dependent transport phenom- 
ena for a quantum dot connected to ferromagnetic leads 
IMHyiHB 12, 13, 14, 15, 16, 17J, giving interesting exam- 
ples in a rapidly developing field of magneto-electronics 
or spin-electronicsdl Hi llfl IM El • 

When several quantum dots are connected to each 
other, they provide remarkable phenomena due to the in- 
terplay of electron correlations, interference effects, etc, 
which depend on how the dots are arranged: e^ double 
quantum dots coupled in series or parallel 23, 24, [2^ |2a, 
iL 28. 22, .30.], T-shaped double dots^^^], etc. In 
particular, the Kondo effect in double quan tum dot sys- 
tems has been observed experimentally i29| i3u], opening 
prospects for exploring further exotic phenomena due to 
electron correlations in nanoscale systems. 

In this paper, we study spin-dependent transport 
through a system of double quantum dots coupled to 
two ferromagnetic (FM) leads in series. Since the Kondo 
effect is sensitive to the presence of the FM leads, a 
variety of remarkable phenomena are expected to ap- 
pear in transport at low temperatures. By applying the 
slave-boson mean-field approximation to the highly cor- 
related Anderson model, we investigate the equilibrium 
and nonequilibrium transport with the use of Keldysh 
Green function techniques. We show that the strength 
of the intcr-dot tunneling is an important key parameter 
to control transport phenomena via the modified Kondo 
resonances. We further discuss how the modified Kondo 
resonances affect transport properties for up- or down- 
spin electrons in the parallel and anti-parallel configura- 



tions of the FM leads. 

The organization of this paper is as follows. In the next 
section, we describe the model and the method. Then 
in §3, we present the results for two different configu- 
rations: spin-polarization of two FM leads is parallel or 
anti-parallel. A brief summary is given in the last section. 



II. MODEL AND METHOD 



We consider spin-dependent transport properties of a 
system with double quantum dots coupled in series. In 
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FIG. 1: Schematic picture of the double dot system coupled 
to the FM leads: tc is the interdot coupling and Vaa{a = 
L, R, a =1, 1) represents the dot-lead tunneling for an electron 
with spin a. 



Fig^we schematically show our double-dot system. It is 
assumed that the left and right leads are spin-polarized, 
where the parallel (P) and anti-parallel (AP) configura- 
tions are considered. In the following discussions, the in- 
tradot Coulomb interaction U is sufficiently large, so that 
double occupancy at each dot is forbid den, which allows 
us to use slave-boson representation^^ |33 of correlated 
electrons in the dots. In terms of this representation, we 
can model the system in FigQ] with a A^(=2) fold de- 
generate Anderson Hamiltonian including an additional 



term for interdot coupling, 



where Fq is the value a,t p — 0. 



H 



/ J ^kacr'^k^a^k^cr + / ^ ^acrJaaJao 



ka ,0" 



TV 



4Y.(fl^''LbyB.. + fLbBblha) 



nrf / J ^a<7\^kaCr^OiJa(7 + Jaa'^a^kacr) 



ka -O" 



E^"(E/-/-+&i^^a-i), 



(1) 



k ay^ka^tr) is the creation (annihilation) operator 



where cj, 

for an electron with spin a in the lead a {a = L, R). Ac- 
cording to the slave-boson representation, the creation 
(annihilation) operator of electrons in the left or right 
dot, d'l^idacr), is replaced by d'l^ -^ fla^a {daa -^ 
faa^a), whcrc 5q {faa) is the slave-boson (pseudo- 
fermion) annihilation operator for an empty state (singly 
occupied state) . The last term with the Lagrange multi- 
plier Xa is introduced so as to incorporate the constraint 
imposed on the slave particles, J2cr=i,i flu faa + ^a^a = 
1. For simplicity, we consider the symmetric dots with 
El^ — Eji^ ~ Eq. The mixing term Vaa in ^ leads to 
the linewidth function 



F^(e) = ^V|y„^,|2<5(£ -£,„,), 



k 



(2) 



which is reduced to an energy-independent constant F" 
in the wide band limit. Here, we introduce the effec- 
tive spin-polarization stren gth p, following the definition 
given in the hterature: [H 111 Ili El El 
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p 






(0 <p< 1), 



(3) 



This quantity is convenient to represent how large the 
strength of the spin polarization is. From the definition 
of p, F" is expressed in the case of the P and AP config- 
urations. 
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To analyze our model, we use the slave-boson mean- 
field approximation jH, IH IH H HI], in which 
boson fields are approximated by their mean values, 
ba{t)/VN -> {ba{t))/VN = ba. By introducing the 
renormalized quantities Vaa — Vaaba, tc = tcb^bji, and 
£a = £o + Aa, the total Hamiltonian JQ) is now replaced 
by 



H 



/ , ^kao-'^k^a-^kacr "1" / , ^aJacrJaa 
ka,a a,<y 

^c2_^\fLafR<J + iRafLa) 
a 

ka -a 



(6) 



To determine the mean-field values of {bL,bR, Xl, Xr), 
we employ the equation of motion method with the non- 
equilibrium Keldysh Green functions. This gives the fol- 
lowing set of self-consistent equations, 



^L(fl) ^ * 



'E 



—r< (r) - - 

^^'-'L,LiR,R)A^) - 2' 



{£l{r) - £Q)b\{R) 



5E/i^(^-^~i(«))^ 



< 

L,L(R.,R) 



(7) 



.(£)=0, (8) 



where G<^(^^^)^(e) 



is the Fourier transform of the 



Keldysh Green function G< „,^(t-i') ee i{fl,{t')f^,{t)). 
The first equation represents the constraint imposed on 
the slave particles, while the second one is obtained from 
the stationary condition that the boson field is time- 
independent at the mean-field level. From the equation 
of motion of the operator faa |25l l3q | , we have the ex- 
plicit form of the Green function. 



^LX(R,R)a 



ie) 



2* {fLiR.){s)f^'^"^ [{S - Sr^D? + (f?^''^) 



+//?(L)(£)f5''^t;'} 



(e - £L + ^f^){e - sr + ttj^) - 4 
I 



(9) 



with t"^ = blT^ (a ^ L,R). Note here that we have de- 
rived the above formulae by applying the analytic contin- 
uation rules to the equation of motion of the time-ordered 



Green functions along a complex contour. 

By using the renormalized parameters determined self- 
consistently, we obtain the current / through the two 



dots, 123, laa 12313 

2ero 



1 = 



J2jdeih{e)-fn{e))T,ie) 



(10) 



with the Fermi function fai^), where the transmission 
probabihty Ta{s) for spin a is 



r.(e) 



24'f^f« 



(e - £L + ifi)(£ - ER + iTS) - 4 



(11) 



We can thus compute the differential conductance G = 
dl/dV from the current / in I|1U|) . Moreover, we obtain 
the density of states (DOS) for the dot a {a — L or R) 
as 



/Oq,o-(£) 



TT 



(12) 



where G^ ^^.{e) is the Fourier transform of the 
Keldysh Green function G^ ^^{t — t') = —i6{t — 
t'){{Uait), fLit')}), which is given by 



G 



^'^'"'"^^^^ " (s - e, + zf.-)(e - .-. + .f«) - t? -^''^ 



This completes our formulation of the non-equilibrium 
transport for the double-dot system with spin polariza- 
tion. 



III. NUMERICAL RESULTS 

In this section, we discuss the results computed at ab- 
solute zero for the P and AP configurations separately. 
For simplicity, the band width of the leads is taken as 
D — GOFq and the bare level for two dots is fixed as 
Eq = — 3.5Fo (Kondo regime). We will use Fq as the unit 
of the energy and choose the Fermi level at ^ = as the 
origin of the energy. 




FIG. 2: Linear conductance Gv=o = dl /dV\v^o as a function 
of the spin-polarization strength p for tc = 0.3,0.5, 1.0 and 
1.5. 
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FIG. 3: The contribution of up and down spin electrons to 
the conductance for (a) tc = 0.3 and (b) 1.5. 



A. Parallel configuration 

Let us start our discussions with the parallel (P) con- 
figuration of the polarized leads. 



1. Equilibrium case [V — 0) 

We first investigate transport properties of the equi- 
librium dot system where the bias voltage V is infinites- 
imally small. In Fig. [3 we plot the linear conductance 
Gv=o — dl /dV\v=a a-s a function of the spin polarization 
strength p for various interdot couplings. For tc = 0.3, 
Gv=o exhibits a maximum around p = 0.7, whereas it de- 
creases monotonically for f c = 1-0 and 1.5. To clearly ob- 
serve the characteristic p-dependence in the conductance. 



let us consider the contribution of up- and down-spin 
electrons to Gv=o- For small interdot coupling {tc — 0.3) 
in Fig. OJa), it is seen that down-spin (i.e. minority- 
spin) electrons contribute dominantly to Gv=o at large 
p, featuring a maximum structure around p — 0.7. How- 
ever, for larger interdot couphngs (e.g. tc — 1.5) in Fig. 
I^l^b), Gv=o gets dominated by up-spin (i.e. majority- 
spin) electrons as p increases. 

The origin of the above characteristic properties can be 
seen in the formation of the modified Kondo resonances 
in the presence of p. Figures 0] and show the local 
DOS for up- and down-spin electrons of the dots around 
the Fermi energy e = 0. In the equilibrium case with 
y = 0, we denote the DOS for the left and right dots as 

Pcrie) = pLa{e) = PRais). 

In the absence of p, one can see the formation of the 
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FIG. 4: DOS for (a) up- and (b) down-spin electrons of the 
dots in the case of tc. — 0.3. The energy e is normahzed by 
the Kondo temperature in the equihbriurn case T^ (T^ ~ 
lO-'Fo). 



Kondo resonance around e = 0, which is rcnormalized 
about 10~^ smaller than the bare resonance. First we 
consider the case of tc — 0.3. As the bare resonance 
width for up-spin electrons increases with the increase of 
p, according to Eq. @ , the resulting Kondo resonance for 
Pl{e) in Fig. Efa) becomes broad, so that P|(0), which 
determines Gv=Oi decreases. On the other hand, as seen 
from pi{£) in Fig. EJb), the width of the Kondo resonance 
for down spin electrons becomes smaller with increasing 
p. Therefore, the corresponding peak in pi{e) splits into 
two sub-peaks when its width gets comparable to the 
renormalized interdot coupling tc- Namely, the Kondo 
resonances in the left and right dots are combined into 
the bonding and antibonding Kondo resonances, which 
are respectively located around 



£± = 77 <! (£l + Sr) ± \l{eL - SrY + 4t ^ 



(14) 



Since Sl = £r &t V = 0, the energy difference 2tc{= 
£+ — £_) gives the splitting of the Kondo resonances. 
Thus, /9x(0) increases first and then decreases with the 
increase of p. The above p-dependence of the Kondo 
resonances naturally explains the characteristic behavior 
of Gv=o in Fig. EJa). 

We next discuss the case of ic = 1-5. As seen in 
Fig. |S1 the large interdot coupling causes the splitting 
of the Kondo resonances even at p = 0, which is ob- 
scured (sharpened) for up spins (down spins) with the 
increase of p. As a result, /O|(0) little changes, while 
Pl(0) monotonically decreases as p increases. This clar- 
ifies why the contribution of up-spin electrons to Gv=o 




FIG. 5: DOS of the dots for (a) up- and (b) down-spin elec- 
trons in the case of f c = 1.5. 



is almost independent of p whereas that of down spin 
electrons decreases with increasing p in Fig. |3fb). 

When we change the strength of tunneling ic, the sys- 
tem changes gradually between the above two typical 
cases. 



2. Nonequilibrium case [V ^ 0) 

We next show the results in the nonequilibrium case. 
By using Eq. H1Q(I , we calculate the current / through the 
system as a function oiV: we set the chemical potentials 
jiL = V/2 and p,R = —V/2. We wish to mention here that 
Aguado et al. studied transport properties for normal 
leads with finite bias voltage. |25l| 

Figure |S1 shows the I-V characteristics as well as the 
contribution of up- and down-spin electrons. Note that 
the conductance Gv=g discussed above coincides with 
the gradient of the /- V curves at T^ = in Fig. For 
tc — 0.3, as shown in Fig. Efa), the current once increases 
and then decreases after making a peak structure in its 

V dependence. The peak structure becomes sharp with 
the increase of p, which is mainly caused by down-spin 
electrons, as seen in Fig. El^b). We can thus see that for 
finite p the current is dominated by down spin electrons 
for small V ^ while it is mainly contributed by up-spin 
electrons for large V . In contrast, we encounter different 
behavior for tc = 1.5 in Figs. EJc) and (d). In the volt- 
age regime shown in the figure, the current / is almost 
linear in V , and is dominated by up-spin electrons as p 
increases. 

As in the equilibrium case, we interpret the above I- 

V characteristics in terms of the modified Kondo reso- 
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FIG. 6: (a) I-V characteristics and (b) the contribution of 
up- and down-spin electrons for tc = 0.3. Similar plots are 
shown in (c) and (d) for tc — 1.5. 
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(c)t,=1.5 
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FIG. 7: DOS for left and right dots at p = 0.8 in the nonequi- 
libriuni case: tc — 0.3 for (a) up- and (b) down-spin electrons; 
tc = 1.5 for (c) up- and (d) down-spin electrons. 



nances. In Fig. [7|we show the Kondo resonances in the 
nonequihbrium case by exploiting tc — 0.3 and tc — 1.5 
with a typical spin-polarization p = 0.8. 

Let us first look at the Kondo resonance at tc = 0.3 in 
Figs. Cl^a) and (b). At V = 0, the DOS for up-spin elec- 
trons shows a broad hump, while that for down-spin elec- 
trons has sharp peaks with small splitting due to tc- With 
increasing V, the Kondo resonance of the left (right) dot 
is shifted, following change in the chemical potential of 
the left (right) lead. Roughly speaking, the current / de- 
pends on the overlap of these two Kondo resonances in 
the energy region between ^l and hr {—V/2 < e < V/2 
in Fig. [TJ. By this reason, at small V, the current / for 
up-spin electrons gradually increases with the increase of 
V and goes down reflecting the decrease in the overlap of 
the Kondo resonances for larger V. On the other hand. 



for down spins (Fig. Efb)), the splitting of the resonances 
is enlarged as V increases, leading to the well-separated 
peaks. Correspondingly, the current I rises sharply for 
small V, and then goes down quickly as V gets large. 

When the inter-dot coupling is large {tc — 1.5), we can 
see the splitting of the Kondo resonances both for up- 
spin and down-spin electrons at V^ = 0, although they 
are well separated (Fig. El^d)) for down spins but more 
or less obscured for up spins (Fig. Efc)). A remarkable 
point is that their shapes little change with increasing 
V, because the coupling between the two dots is much 
larger than the applied bias. Note that the DOS for up- 
and down-spin electrons is almost flat in the vicinity of 
the Fermi energy, so that the current / is approximately 
linear in the bias voltage and mainly caused by up-spin 
electrons, in accordance with the results of Figs. EIc) and 
(d). Here, we wish to mention what happens for larger 



voltage in the case of tc = 1-5. We have checked that 
the system undergoes a first-order phase transition when 
V further increases beyond V/Tk ''^ 4, where the cur- 
rent drops discontinuously. This imphcs that a strongly 
coupled molecule is formed between two dots, which is 
almost decoupled from two leads. Such a phenomenon 
was already pointed out in the p — case. l25| The detail 
of the first-order transition can be found there. 



B. Antiparallel configuration 

We now move to the system with the AP configura- 
tion of the polarized leads. We can continue discussions 
similarly to the above P case, so that we summarize the 
essential points briefly. 



1. Equilibrium case {V — 0) 
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FIG. 8: Linear conductance Gv^o as a function of the spin 
polarization strength p for the AP configuration. 

We start with the equilibrium case. In Fig. |S1 we 
plot the linear conductance Gv=o ^ a function of the 
spin polarization strength p. Note that the contribution 
of up- and down-spin electrons to the total conductance 
is same in the AP case. Nevertheless, we observe the 
behavior somewhat analogous to the case of the P con- 
figuration: for small tc the total conductance features 
a peak structure in its p dependence, while it decreases 
monotonically for large tc- On the other hand, for the 
fully polarized case with p = 1, the current vanishes in 
the present case in contrast to the P configuration. 

Similarly to the P case, the characteristics in the con- 
ductance can be interpreted in terms of the modified 
Kondo resonances. In Fig. |51we show the local DOS. Ac- 
cording to Eq.©, we see that pl]{£) [pLii^)) is identical 
to PRi{e) (pfl,|(e)). We focus on the up-spin electrons for 
simplicity, pri (e) and pLi (s) are mixed with each other 
in the presence of the coupling tc- We note that the right 
(left) resonance gets narrower (broader) with the increase 
of p, so that for tc = 0.3 there appears a dip structure 
in PL^i^) around e = 0, while a single peak structure 




(d)t,=1.5 ___^0.3 

/■"■"■\ p=0.6 

./.. .A,. p=0.9 




e/r/ 



FIG. 9; DOS of the left and right dots at various p: tc = 0.3 
in (a) and (b); tc = 1.5 in (c) and (d). 



still remains in pfl|(e), as seen from Figs. EJa) and (b). 
Recall that the conductance is roughly given by the prod- 
uct of PR.lis) and pLiis) at the Fermi level. Therefore, 
we naturally see that the combination of the peak and 
dip structures in the DOS gives a maximum structure 
in the conductance for tc = 0.3. On the other hand, for 
tc — 1.5, the development of the peak structure in p/j|(£) 
is slow with the increase of p. Thus the development of 
the dip structure in pLiis) dominates the slight change 
in p/j|(e), leading to the monotonic decrease of the con- 
ductance as a function of p, as seen in Fig. |H1 Note that 
at p = 1, the current vanishes, in contrast to the P con- 



figuration, because either of the DOS of the right or left 
dot vanishes at the Fermi level. 



2. Nonequilibrium systems (V / 0) 

We finally show the results in the nonequilibrium case 
with the AP configuration. Figure [TUl shows the I-V 
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FIG. 10: /- y characteristics for the AP configuration at p = 
0,0.4 and 0.8: {a)tc = 0.3 and (b)fc = 1.5. 

characteristics of tc — 0.3 and 1.5 for various choices of 
p. We note again that the ^-dependence of the total 
current shows similar behavior to that for the P configu- 
ration (Fig. O although the origin of the characteristics 
is somewhat different, as mentioned above. Shown in 
Fig. ^3 are the Kondo resonances in nonequilibrium case 
at p = 0.6. When the bias voltage is applied, the peak 
position of pi,o-(e), PRa{£) (ct = t or |) follows the change 
of the left or right chemical potential, as shown in Figs. 
Illf a) and (b). As mentioned above, the current for up- 
or down-spin electrons is controlled by the overlap of the 
Kondo resonances of the left and right dots in the region 
of e = -F/2 to £ = T//2 in Fig. El This gives rise 
to the maximum structure in the /- V characteristics for 
tc — 0.3. As seen in Figs. Illf c') and (d) for tc = 1.5, 
the shape of the Kondo resonances with splitting is not 
altered so much up to a certain voltage, giving a linear 
conductance, as is the case for the P configuration. For 
larger voltage, the system suffers from a sudden phase 
transition accompanied by discontinuous drop in / |25| . 
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FIG. 11: DOS of the left dot (solid line) and right dot (broken 
line) for the bias voltage V — and 3T^ at p = 0.6: (a) up- 
and (b) down-spin electrons in the case of tc = 0.3; (c) up- 
and (d) down-spin electrons in the case of tc = 1.5. 



IV. SUMMARY 

We have discussed equilibrium and non-equilibrium 
transport properties of double quantum dots attached to 
the FM leads in series. It has been found that, under the 
influence of the FM leads, transport shows the charac- 
teristic properties due to the modified Kondo resonances. 
We have clarified that the total conductance shows anal- 
ogous p-dependence between P and AP configurations of 
the leads, although the origin of such behavior is some- 
what different from each other. 

For the P configuration, at small inter-dot coupling 
tc, the conductance of down-spin (minority spin) elec- 
trons dominates that of up-spin (majority spin) elec- 
trons, giving rise to a peak structure in the p-dependent 
conductance. On the other hand, for large coupling tc, 



the conductance exhibits monotonic p-dependence, and 
is mainly contributed by up-spin electrons in the large p 
region. For the AP configuration, we have observed sim- 
ilar p-dependence in the total conductance, although the 
contribution of the up- or down-spin electrons is always 
same in this case. 

We have also studied the current / at finite bias volt- 
age, which again shows distinct properties depending on 
tc- For small inter-dot coupling, the Kondo resonances of 
the left and right dots are shifted, following the change 
in the chemical potential. This gives rise to a charac- 
teristic peak structure in the V dependence of the cur- 
rent. On the other hand, for large inter-dot coupling, the 
shape of the modified Kondo resonances is not changed 
so much up to large V, leading to the linear conductance 
in the wide V region, beyond which the systems under- 
goes a first-order transition accompanied by the sudden 



decrease of the current. 

Experimentally, the Kondo effect in double dot systems 
has been observed recently ^9, 30). We thus expect that 
in the near future, it may become possible to observe the 
Kondo effect in double dot systems coupled to the FM 
leads, providing further interesting examples of correla- 
tion effects in the context of spin-dependent transport in 
nanoscale systems. 
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